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For a dynamical map A(t, 0), which sends a state p(0) of quantum open system to a state p(t) = 
A(i,0)p(0), the decomposition law A(£, 0) = A(£, £ C )A(£ C ,0) may break down at a specific time t c . 
In this paper, we present a method to find the singular points t c and propose a measure for the 
singularity of the dynamical map. Two examples are portrayed to illustrate the method, the measure 
of singularity for these singular points is calculated and discussed. An extension to high-dimensional 
system is presented. 
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I. INTRODUCTION 

The actual dynamics of any real open quantum system 
is expected to deviate to some extent from the Marko- 
vian evolution. This deviation can be measured by non- 
Markovianity and it has attracted much attention in re- 
cent years, leading to a deeper understanding of quite a 
few issues in the theory of open quantum system [l|-|8( . 

Non-Markovian systems can be found in many 
branches of physics, including quantum optic s M , Il0| . 
solid state physics [ll|, quantum chemistry pjj], and 
quantum information processing [13| . Since non- 
Markovian dynamics modifies monotonic decay of quan- 
tum coherence, it may protect quantum entanglement 
in composite systems for longer time than standard 
Markovian evolution 1A\. In particular it may protect 



the system against the sudden death of entanglement 
|15| . Therefore, it is interesting to quantify the non- 
Markovianity within the description of quantum open 
system. 

There are two approaches to quantify the measure of 
the degree of non-Markovianity. One approach is based 
on the idea of the composition law which is essentially 
equivalent to the idea of divisibility [16j . This approach 
was used recently in Ref . [171, llll to construct the mea- 
sure of non-Markovianity, quantifying actually the devi- 
ation of the dynamical map from divisibility. Another 
approach is as in Ref. [19j | , where the authors define non- 
Markovianity dynamics as the information flow from the 
environment back into the system, the measure manifests 
itself as an increase in the distinguishability of pairs of 
evolving quantum states, and the information is identi- 
fied to be the Fisher information [20| . 

The measure for non-Markovianity proposed in 
Ref.[17| is based on the completely positive divisibility 
of a dynamical map: a trace preserving completely pos- 
itive map A(i,0) is completely positive divisible (CP- 
divisibility) if it can be written as, 



and A(i2,ii) is completely positive for any ti and t\ 
(£2 > tj. > 0). By contrast, we say that the map A(t2,0) 
is positively divisible (P-divisibility) if A(i 2 ,£i) sends 
states into states but it is only positive, and that A (£2, 0) 
is indivisible if neither P-divisibility nor CP-divisibility 
holds. 

In this paper, we shall consider the other situation 
where, 



A(t 2 ,0)?6A(t2,*c)A(tc,0), 



(2) 



at a special time t c , £2 < t c < 0. We will refer to this 
instance of time t c as the singular point of the dynam- 
ical map A(£ 2 ,0),t2 G (0,oo). Taking a qubit (two-level 
system) as an example, a method to find the singular 
point is presented, a measure to quantify this singularity 
is proposed and discussed. 

This paper is organized as follows. In Sec. II, we present 
a general formalism for a qubit dynamics, exhibiting the 
method to find the singular point t c . A measure to quan- 
tify the singularity is constructed. Two examples, one 
describes a qubit coupled to a harmonic oscillator bath 
and the other includes a qubit coupled to a finite spin 
bath, are given to illustrate the critical point in Sec. Ill, 
the measure of singularity is also calculated and discussed 
in this section. A generalization of the representation to 
(i-dimensional open systems is presented in Sec. IV. Fi- 
nally, we conclude our results in Sec. V. 



II. GENERAL FORMALISM FOR A QUBIT 
DYNAMICAL MAP 



Consider a dynamical map A(£, 0) for a qubit (or 
two- level system), which sends an arbitrary initial state 
p(0) = (1 + n(0) ■ ct)/2 with Bloch vector n(0) = 
(n x (0),n y (0),n z (0)) into a state p(t), 



A(i 2 ,0)=A(i 2 ,ii)A(ii,0), 



(1) 



p(t)=A(t,0)p(0)=-(l + n(t)-a). 



(3) 



Without loss of generality, the Bloch vector n(t) — 
(n x (t),n y (t),n z (t)) can be written as, 

n(t) = n(0) ■ D(t) + f(t), (4) 

where D(t) is a 3 x 3 matrix and f(t) is a time-dependent 
vector. 

Now we elicit the condition for A(t 2 ,0) 7^ 
A(t 2 ,t c )A(t c , 0). To this aim we introduce an ancilla A 
and define, 



and 



M s a = A(t 2 ,*c) <8> Ia(|$sa}($sa|), 



(5) 



where |$sa) = |0)s®|0)a+|1)s®|1)a is an unnormalized 
maximally entangled state of the qubit and ancilla, and 
Ia denotes the identity operator of the ancilla. Since the 
ancilla is also a qubit, Ms a can be written as, 



M SA = -{x + r-a^ + s-as + Ss-V -Ua) 



(6) 



Here a; is a constant, r and s are vectors, I is an identity 
x f 



matrix, F = I "Z, J . 1 and V is a 3 x 3 matrix, which 

is determined by the map A(t2,t c ) an d will be derived 
in the following. If A(i 2 ,0) = A(t 2 ,t c )A(t c , 0) holds, the 
map A(t 2 ,t c ) would send the state p{t c ) = A(t c , 0)p(0) to 
state p(t 2 ). In terms of Msa, this can be expressed as, 



p(t 2 ) = A(t 2 ,t c ) P {t c )^Tv A [MsAls®p T A {tc)] 



(7) 



s(t 2 ,t c ) = f(t 2 )-f(t c )-V. 



(12) 



From the above derivations, we find that once x, f, 
s, V are (uniquely or non-uniquely) established for 
any t c {t 2 > t c > 0), the decomposition A(i 2 ,0) = 
A(t 2 ,t c )A(t c ,Q) holds true, namely there are no singular 
points in the time interval [0, t 2 ]. We notice that the null- 
determinant of D(t c ) plays a key role in finding V, it can 
thus be taken as a condition to find the singular point t c 
if the matrix D(t 2 ) is of full rank. Mathematically, the 
necessary and sufficient condition for Eq. (|10[) to have 
no solution is that the rank of D(t c ) must be smaller 
than the rank of [D(t c )\D(t 2 )}, where [D(t c )\D(t 2 )} is an 
augmented matrix obtained by attaching the columns of 
D(t 2 ) to the columns of D(t c ). 

To quantify the singularity of the singular point, we 
introduce the trace distance, D(pi,p 2 ) = |Tr|pi — p 2 \, 
which is an appropriate measure for the distinguisha- 
bility between two quantum states pi and p 2 . Here 
\A\ = VA^A. We define the singularity measure of a dy- 
namical map A (£2, 0) at time t c by 

S A (t c )=maxD(p(T),p te (T)), (13) 

where p(T) = A(T,0)p(0) and p tc (T) = A(T,t c ) ■ 
A(t c , 0)p(0) is the solution of the dynamical process tak- 
ing p(t c ) = A(t c ,0)p(0) as initial state. The maximum is 
taken over all initial states and the final time T. 



III. EXAMPLES 



Writing p(t 2 ) = |(I, as) ( -./. \ ) , we obtain from 
Eq.©, 



1 
n{t 2 ) 



F 



n(t c ) 



x + f ■ n(t c ) 
s + H(t c ) ■ V 



(8) 



It is easy to find that, 



X = 1, r = 0, 
n(t 2 ) = n(t c ) ■ V + s. 



(9) 



Considering that p(0) is an arbitrary initial state, namely 
n(0) is arbitrary, Eqs. (|4|9|l together yield, 



D{t 2 ) = D(t c ) ■ V, 
f(t 2 ) = /(t c ) -V + s. 



(10) 



The condition for A(£ 2 ,0) 7^ A(t 2 ,t c )A(i c , 0) now is 
equivalent to that there does not exist a matrix V to 
satisfy Eq. (fT0|) . If the determinant of D(t c ) is non-zero, 
we have 



V = D- 1 (t c )-D(t 2 ), 



(11) 



In this section, we will present two examples to illus- 
trate the singular point and the measure of singularity. 
The first example is a dephasing model that consists of 
a spin-i particle coupling to a spin-bath. The coupling 
Hamiltonian commutes with the free Hamiltonian of the 
central spin, thus the central spin conserves its energy. 
In the second example, we consider a dissipative system, 
the energy of the system is no longer conserved. 



A. A two-level system coupling to a finite spin bath 

Consider a central spin- i coupling to a bath of N spin- 
\ particles. The interaction Hamiltonian is, 



H 



N 

J2 A ^z<J k zl 
fc=i 



(14) 



where A^ — A/y/N represents the coupling constants. 
Assume the initial state of the whole system is p s (0) <§) 
(w-Oi i- e -> a H spins in the reservoir are in a maximal 
mixed state. The density matrix of the central spin at 



time t takes, 



p(t) 



PXX Pl2C0S-(if) 



P 2 i-s^(^) 



P22 



(15) 



In terms of dynamical map, the dynamics can be rep- 
resented as, A(t,0)p = |(1 - cos N (^M))a z pa z + |(1 + 



.,N (2At 

S ^ 

equation, 



cos ( :f M=))p. This is equivalent to the following master 



p = ~f(t)£(p), 



(16) 



where C(p) = cr z pa z — p, and the time-dependent decay 
rate is 7(f) = A-\/]Vtan(^4=). This model is discussed 
in several papers as a typical example to quantify non- 
Markovianity. 

Writing pit) in Eq. (|T5|) in the form of Eq. (|4]), we find 

C(t) ' 
D(t) = ( C(t) I , 
1 



/(*) = 0, 



(17) 



where C(t) — cos JV (^^). The singular point t c can be 
found by solving C(t c ) — 0, it yields, 



N 

L = lZ (2n + 1)7r ' n = ' 1 ' 2 '"- ( 18 ) 

By the definition of the measure of singularity, we obtain, 
S A (t c )= max |C(T)||p 12 (0)|, (19) 

|P12|,T 



of the initial density matrix p(0) = 



where T is a time, T > t c , and pi2(0) denotes the element 

Pn(0) pi 2 (0)' 

P2l(0) p 22 (0) 

After a simple algebra, we have Sa(£ c ) = \ for any sin- 
gular point given in Eq.(fT5|). It is interesting that the 
singularity measure of these singular points are equal. 
Indeed, examining the dynamical map A(£, 0), we find 
that the features of A(t, 0) around any t c are the same. 



the creation and annihilation operators of k — th oscilla- 
tor, which couples to the system with coupling constant 

9k- 

This model is exactly solvable [9j. Assuming the sys- 
tem and the reservoir initially uncorrelated, we can ob- 
tain a time-dependent master equation in the interaction 
picture, 



At), , 

1 



j(t)(a pA 



-a' a p- -po a 



(21) 



where e(t) = -2Im[|£$] and j(t) = -2Re[|^]. e(t) 
plays the role of Lamb shift and j(t) is the decay rate. 
Both e(t) and j(t) are time-dependent. c(t) is deter- 
mined by c(t) = — L f(t — T)c{T)d(r) 1 where f(t — r) = 
J dcdJ(uj)exp(i(ujQ — uj)(t — T)) is the environmental corre- 
lation function. In the derivation of the master equation, 
the reservoir is assumed in its vacuum at t = 0. 

Consider the following spectral density, J(ui) = 
— / 2A + x 2 ' wnere 7o represents the coupling constant 
between the system and reservoir, A defines the spec- 
tral width of the coupling at the resonance point u)q. 
For the spectral density J(w), we have e(t) = 0, c(t) = 
c e~ 



- A */ 2 [cosh(f ) + |sinh(f )], and 



lit) 



27 Asinh(dt/2) 



dcosh(di/2) + Asinh(di/2) 



(22) 



initially in p(Q) = 
Hamiltonian i 
at time t, p(t) 



, by the effective 



with d = •*/ A 2 — 270 A in Eq. ([2~Tj) . Assume the system 

Pee(0) Peg(0) 
PaeiO) P 3 M 

Hamiltonian approach [2 lj. we have the density matrix 

Pee{t) Peg(t)\ ^^ 
Pge(t) Pgg(t) J ' 

Pee(t) = p e Mz~ ! °~< {t ' )dt \ 
Pgg(t) = 1- p ee (t), 

Peg(t) = p; e (t)=e-if«^ dt 'p eg (0)- (23) 



B. The damping J-C model 

This example consists of a two-level system coupling 
to a reservoir at zero temperature. The reservoir con- 
sists of infinite number of harmonic oscillators that is 
also referred in the literature as the spin-boson model. 
The Hamiltonian for such a system reads, 



H = H + H T , 



(20) 



where H = fkOoa+a^.+Ylk^kbkbk, Hi = (j + B + (j-B\ 
and B = J2k9k^k- The Rabi frequency of the two- level 
system and the frequency for the k — th harmonic oscilla- 
tor are denoted by wq and W&, respectively. b\ and bk are 



It is easy to show that the matrix D{t) and fit) in this 
example are, 

( D u (t) \ 

D(t) = DaaC*) (24) 

V D 33 (t)J 

and 

/(*) = (/x, fyJz) = (QA(e- f ^ (t ' )dt ' - 1)), (25) 
where 



Dn = D 2 2 = e 2 
#33 = D\ x . 



1 /o 7(*')*' 



(26) 



We find from Eq. ([24]) that D j:j (t c ) = 0, j = 1,2, 
or 3, gives the singular points. Djj(t c ) = can hap- 
pen only when 70/A > 1/2. Noticing that Du(t) — 
e- At / 2 [cos(^) + ^sin(^p)] for 7o /A > 1/2, where 

do = \/\X 2 — 2toA|, we obtain the nth singular point 



t. 



(n) 



= iM-H- 



.) + n7r),n = 0,1,2, 



At 



-y/|l-2 7o /A| 

these singular points, the singularity measure can be 
given by maximizing the distance D(p(T), p tc (T)) = 
yD 33 (T)(n((0) +nf(0)) +D* 3 (T)(1 + n 3 (0)) 2 over T 
and the Bloch vector n(0) = (ni(0),n 2 (0),n 3 (0)) with 
constraint < n 2 (0) + n 2 (0) + n|(0) < 1. Simple alge- 
bra shows that the maximum of the nth singular point 



arrives at T 



rp(n) 



2(w+l)7T 

do 



0,1,2..., and 



«3(0) = T^gS),n 2 (0) + nl(0) = l-n 2 (0). The mea- 

(n) 

sure of singularity for the nth singular point t c is then 



5a (* 



(n)i 



A 



ny-r, for < D 33 (T^) < 0.5. 



SaW) 



(n)\ _ 



1-e" 

2(i + 1)ttA 



d 



, for 0.5 < D 33 (r (n) ) < 1. (27) 



We find that the values of singularity measures are dif- 
ferent in the two examples. In the first example, the sin- 
gularity for all singular points are the same, Sa(£c) = \, 
while in the second one, the singularity depends on the 
singular points. This results from the difference in the 
states at the singular point t c in the two examples. Es- 
pecially, as n increases, the singularity decreases and fi- 
nally tends to zero as n — > 00. In other words, the sin- 
gularity of larger t c is smaller than that for a smaller 
t c . This can be understood as that at large t c , the state 
of the open system is more close to the steady state, 
leading to a small difference in the states. Furthermore, 
the difference in singularity is a reflection of the system- 
environment coupling. The first example is a dephasing 
model, it conserves the system energy but spoils the off- 
diagonal elements of the density matrix. By contrast, the 
system would decay to its ground state at the singular 
points in the second example. 

Since the measure proposed here quanti- 
fies the non-divisibility of the dynamical map, 
A(t, 0) 7^ A(t,t c )A(t c ,0), hence it can measure the 
non-Markovianity of the map. The non-Markovianity 
in this situation depends both on the measure of the 
singularity and the number of singular points. Therefore, 
we propose, 



N u =Y,SdtP), 



(28) 



to quantify the non-Markovianity caused by the singular 
points tc', (j = 1,2,3, ...). Physically, once a dynamical 
map has a singular point t c , the state at time t > t c 
would depend on the state at an earlier time t' < t Cl 
although the state at t c is the same. This feature can be 
found by examining Eq. (|15p . which is a reminiscence of 
the classical non-Markovian process. 



The present prediction can be observed in the exper- 
imental setup in [22j, where the polarization degree of 
freedom of photons plays the role of open system, the 
environment was simulated by the frequency degree of 
freedom with two central frequencies at lj\ and w%. The 
evolution of the off-diagonal elements of the photon den- 
sity matrix takes, \H)(V\ -> K*(t)\H){V\, \V)(H\ -> 
K,(i)\V){H\. Here n(t) is adjustable and can be manipu- 
lated to zero at times t c — — A ^. A „ , n = 0, 1, 2, ..., where 
Alo = LU2 — Wi, is the difference in the central frequen- 
cies of the environment. An denotes the difference in the 
refraction indices of horizontally and vertically polarized 
photons. The observed final states are different that de- 
pend on whether an observation is made at the singular 
points t c . By measuring the difference in the final states, 
the singularity can be quantified in the experiment. 



IV. EXTENSION TO d-DIMENSIONAL 

SYSTEMS 

We consider now an arbitrary dynamical map A(i, 0) 
with t > for a quantum d-dimensional system (qudit). 
Let {A M }™ =1 with n = d 2 — 1 be a set of traceless qudit 

Together with 



dS 



ft u ■ 



observable satisfying Tr^X^Xu) 
the identity operator they form an orthonormal basis 
for all the qudit operators. Thus we have expansions 
g = (I + n ■ X)/d for the initial state and A(t, 0)g = 
(I + n(t) ■ X)/d with n(t) = D(t) ■ n + e(t) for the final 
state g(t) = A(t,0)g. Here ft = Tr(gX), n(t) = Tr(g(i)A), 
and e(t) = Tr[(A(i,0)/)A] are n-dimensional real vec- 
tors and D(t) is an n x n real matrix with matrix ele- 
ments given by [[£>(*)]W = Tr((A(i,0)A„)A M ). The lin- 
ear trace-preserving map A(t, 0) is determined uniquely 
by D(t) and e(t) and vice versa. For later use we denote 
by V(t) = {ae R n \D{t) ■ a = 0} the null space of D(t). 
Let t > t c > and consider the possible decompo- 
sition of a dynamical map A(i, 0) = A(t,t c )A(t c ,0) for 
some linear trace-preserving map A(t,t c ). Any linear 
qudit map, e.g., A(t,t c ), is in a one-to-one correspon- 
dence with a 2-qudit operator, e.g., R = A(t,t c ) <S>Z(<&) 
with $ being the projector of the subnormalized 2-qudit 
state |$) = J2 n \ n ' n )- We denote by S the n x n 
matrix with elements [[S]] M „ = Tr(i?A M <E) X^)/d for 
u, v = 1,2,..., n = d 2 - 1 and f= Tr(i?A® I)/d. For a 
trace-preserving map it holds Tr (RI(& A) = and there- 
fore S and r determine uniquely R and consequently the 
linear trace-preserving map A(i, t c ). By definition the lin- 
ear map A(i, t c ) is a possible decomposition if and only 
if for an arbitrary initial state g it holds 

g{t) = A(t, 0)g = A(t, t c )(A(t c , 0)g) = A(t, t c )g{t c ). 

, (29) 
Lemma. If a qudit operator O satisfies (ip\0\ip) = 

for an arbitrary pure qudit state \ip) then O = 0. 

Proof. Let V12 = Y^ij \hj)(J^\ be the swapping oper- 
ator of 2 qudits and I12 be the identity operator. From 



the following identity 
W 12 r- 



^VWlglVWH ^fi" (30) 



it follows that = d{d+ l)Tri ((Cq ® J 2 )Wi 2 ) = 0. 

Theorem. Given a qudit channel A(i, 0) with t > t c > 
there exists a linear trace-preserving map A(t,t c ) such 
that A(i, 0) = A(t, * c )A(t C) 0) if and only if V(t c ) C V(t). 
Moreover the decomposition A(i, t c ) is unique if and only 
if det£>(i c ) ^ 0. 

Proof. Necessity (only if part), i.e., A(i, 0) = 
A(t,t c )A(i c ,0) infers V(t c ) C V(t). From Eq.® 
it follows that for arbitrary g it holds Tr (Xg(i)) — 
Tr(XA(t,t c )g(t c ))) = Tr(RX ® g T (t c )) = S ■ n{t c ) + f. 

Taking into account ft(t) — Tr(Ag(£)) = D(t) -n + e(t) for 
t and t c we see that (D{t)-S-D(t c ))-n = S-e(t c )+r-e(t) 
must hold for arbitrary ft = Tr(pA) with g being a den- 
sity matrix. If we let ft = with corresponding state 
being g = I/d then we obtain f = e(t) — 5 ■ e(i c ). As a 
result A • ft = 0, where A := D(t) - S ■ D(t c ), for arbi- 
trary ft = Ti(gX) with g being a density matrix. From 
A • Tr(gX) = it follows that Tr(gA, J ,) = for arbitrary 
/i and arbitrary qudit state g where A^ = J2 V A^Aj, is 
traceless, i.e., TrA p =0. As a result of lemma, A M = 
for all fi, i.e., D(t) = S-D(t c ). Thus we have V(t c ) C 7(t) 
since D(t c ) ■ ft = infers D(t) ■ ft = S ■ D(t c ) ■ ft = 0. 

Sufficiency (if part), i.e., V(t c ) C V(t) infers A(*,0) = 
A(t,i c )A(i c ,0). Let {ej^ span V(t c ) where K = 
dimV r (i c ) and {ei}™ =K+1 span the orthogonal comple- 
ment of V(t c ). As a result D{t c ) ■ S% = and thus 
D(t) ■ &i = for i = 1,2,..., if since we have as- 
sumed V(t c ) C V(t). The equation £>(t) = 5 • 
D(t c ) is equivalent to D(t) ■ e, = S ■ D(t c ) ■ Si for 
i = 1, ...n. Since D(t c ) ■ ei = I?(i) • ei = for 
i = 1, ...,K the equation becomes -D'(i) = 5 • D'(t c ) 
with D'(t c ) = [L>(i c )eft: + i,-D(tc)eK + 2,...-D(t c )e„] and 
D'(i) = [£>(t)? x +i,-D(*)gjf+2,-..,I>(*)^,] being of di- 
mension n x (n — if). Since the rank of D(t c ) is n — K 
there are exactly n — K (among n) linearly independent 
row vectors of D'(t c ). Therefore it is always possible to 



expand each row vector of D'(i), an [n — K )-dimensional 
vector, by those n row vectors of D'(t c ), i.e., for any given 
i there exist real numbers Sik such that 

([D'(t)] il ,[D'(t)] i2 ,...,[D'(t)] i , r ^ K ) = 

n 

J2 Sik([D'(t c )]ki, [D'(t c )] k2 , . . • , [D / (t c )] fclTl _jf)<31) 
fe=i 

The n x n matrix S formed by the coefficients Sik with 
i, k = 1,2, ... ,n satisfies D(i) = S ■ D(t c ) and together 
with f = e(t) — S ■ e(t c ) determines R and thus A(t,t c ) 
such that A(t,0) = A(t,i c )A(i c ,0). Moreover if and only 
if K = 0, i.e., V(t c ) is empty, i.e, det D(t c ) ^ 0, S as well 
as r = e(t) — S ■ e(t c ) is unique and therefore A(t,t e ) is 
unique. 

We note that the condition V(t c ) C V(t) is 
equivalent to Rank[£)(i c )|£)(t)] = Rank(Z)(t c )), where 
[D(t c )\D(t)) is an augmented matrix. The decomposi- 
tion A(t, 0) = A(t, t c )A(t c , 0) does not exist if and only if 
Rank[D(i c )|D(t)] > Rank(D(t c )). 



V. CONCLUSION 

In summary, we have explored the singular point t c 
where the dynamical map A(t, 0) 7^ A(t,t c )A(t c ,0), i.e., 
the dynamical map is indivisible at the instance of time 
t c . We quantify the singularity of the singular point t c 
and present examples to show the singularity. Until now 
these points were not aware in the divisibility-based mea- 
sure of non-Markovianity, hence it would contribute to 
the understanding of quantum non-Markovian process. 
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